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fluxes.  These  results  compared  satisfactorily  with  climatological  studies  of  the 
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1 .  INTRODUCTION 


Determination  of  the  turbulent  (latent  plus  sensible) 
heat  flux  between  the  ocean  and  atmosphere  is  crucial  for 
understanding  several  aspects  of  the  earth's  atmosphere- 
ocean  system.  The  turbulent  heat  fluxes  combined  with  the 
radiative  fluxes  determine  the  net  exchange  of  heat  across 
the  air-sea  interface.  It  is  estimated  that  deep  water 
formation  in  the  Greenland,  Norwegian  and  Barents  Seas 
comprises  about  60%  of  the  total  volume  of  the  Arctic  Ocean 
deep  water  (Swift,  Takahashi  and  Livingston,  1983). 
Therefore,  accurate  measurements  of  the  turbulent  heat  flux 
in  this  region  are  necessary  to  understand  deep  oceanic 
convection  and  deep  water  formation. 

The  measurement  of  heat  flux  between  the  ocean  and 
atmosphere  is  also  an  important  factor  in  calculating  global 
circulation  of  the  ocean  and  atmosphere  and  therefore  is 
critical  to  a  complete  understanding  of  the  global  climate. 
Heat  flux  is  directly  linked  to  oceanic  convection  and  deep 
water  formation  in  the  Arctic.  Formation  of  Arctic  deep 
water  causes  southerly  flow  at  depth,  and  the  water  that 
sinks  is  replaced  at  the  surface  by  the  northerly  flowing 
warm  surface  waters  of  the  North  Atlantic.  This  conveyor 
belt  of  oceanic  flow  is  prominent  in  regulating  the  global 
climate.  The  warming  of  the  ocean  surface  by  the  surface 
currents  increases  atmospheric  instability  leading  to  larger 
heat  fluxes  from  the  ocean  to  the  atmosphere  and  increased 
cyclogenesis.  The  increase  in  storms  causes  more 
precipitation,  lower  surface  salinity,  increased  ice 
formation  and  lower  air  temperature.  As  the  sea  ice  melts 
the  surface  water  becomes  more  buoyant,  and  the  water  column 
is  strongly  stratified.  The  stratification  in  turn  hinders 
convective  overturning,  slowing  the  formation  of  deep  water 
and  completing  a  negative  feedback  loop  with  a  cycle  of 
about  20  years  (Bourke,  1994) . 
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Knowledge  of  the  amount  and  variability  in  time  and 
space  of  heat  flux  is  also  needed  to  predict  the  depth  and 
other  properties  of  the  atmospheric  boundary  layer  and  the 
oceanic  mixed  layer.  Turbulent  heat  flux  and  many 
properties  of  the  atmospheric  boundary  layer  and  the  oceanic 
mixed  layer  are  used  in  synoptic  and  mesoscale  atmosphere 
and  ocean  models.  Therefore,  accurate  measurement  of  the 
fluxes  offers  a  chance  for  improvement  in  the  models. 
Additionally,  an  accurate  account  of  heat  flux  is  required 
to  model  the  formation  and  dissipation  of  sea  ice  properly. 

Unfortunately,  the  Arctic  has  relatively  few 
observations  and  is  therefore  less  well  understood  than  the 
mid-latitudes.  Few  measurements  of  sensible  and  latent  heat 
fluxes  have  been  made  over  the  polar  latitudes.  Vowinckel 
and  Taylor  (1964)  (hereafter  VT) ,  Bunker  (1976)  and  Gorshkov 
(1989)  have  made  climotological  studies  of  the  turbulent 
heat  fluxes.  Additionally,  Hakkinen  and  Cavalieri  (1989) 
(hereafter  HC)  examined  turbulent  heat  fliixes  using  gridded 
Fleet  Numerical  Meteorology  and  Oceanography  Center  (FNMOC) , 
Naval  Operational  Global  Atmospheric  Prediction  System 
(NOGAPS)  surface  analysis  for  the  year  1979.  Their  results 
9&ve  extremely  high  values  of  total  turbulent  heat  flux  near 
the  ice  edge  in  the  winter  months . 

In  this  study,  measurements  of  the  turbulent  heat 
in  the  Greenland,  Norwegian  and  Barents  Seas  were 
based  on  direct  measurements  of  bulk  variables.  Usually, 
computed  fluxes  were  within  the  calculated  error  of  the 
climatological  studies.  Additionally,  comparisons  were  made 
with  the  results  of  calculating  the  sensible  and  latent  heat 
flux  based  on  10  minute  averaged  bulk  parameters  versus  ship 
cruise  averaged  bulk  parameters.  These  results  show  the 
^otal  turbulent  heat  fluxes  are  within  5%  irrespective  of 
averaging. 
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In  Chapter  II,  previous  studies  in  the  measurement  of 
heat  flux  in  the  Greenland,  Norwegian  and  Barents  Seas  are 
reviewed.  In  Chapter  III,  the  data  sets  and  methods  of 
collection  are  described.  Chapter  IV  contains  a  description 
of  the  bulk  formula  used  to  compute  the  turbulent  heat 
fluxes.  The  results  of  this  exercise  and  comparisons  with 
previous  studies  are  included  in  Chapter  V  and  conclusions 
are  contained  in  Chapter  VI. 
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II.  PREVIOUS  STUDIES 


In  a  review  of  the  literature  regarding  heat  flux,  four 
sources  were  valuable.  The  first  three  studies  are 
climatological  works,  whereas  the  final  study  emphasizes  the 
relationship  of  heat  fluxes  to  numerical  ocean  and 
atmosphere  modeling.  For  each  study  the  author's  methods, 
data  sources,  and  assumptions  made  in  the  interpretations 
will  be  reviewed. 

A.  VOWINCKEL  AND  TAYLOR 

Vowinckel  and  Taylor  (1964)  calculated  long  wave 
radiation,  latent  (evaporation)  and  sensible  heat  flux  over 
the  Greenland,  Norwegian  and  Barents  Seas  north  of  65°  N. 

Their  calculations  are  after  Sverdrup  (1951) 

E=k(e^-e^)  V  (2.1) 


LHF=E  (2.2) 


where  E  (mm/24  hrs)  is  the  evaporation  rate;  k  is  a 
dimensionless  coefficient  called  the  evaporation  factor 
which  is  empirically  determined  from  an  average  for  the 
entire  area.  VT  used  0.104  for  k;  The  parameters  e,  (mb) 
and  eg  (mb)  are  the  water  vapor  pressures  at  the  sea  surface 
and  at  the  shipboard  measurement  height,  respectively;  V 
(m/s)  is  the  wind  velocity  as  measured  onboard  ship;  LHF 
(cal  /  cm^  month)  is  the  latent  heat  flux;  and  (cal)  is 
the  calculated  latent  heat  of  vaporization  at  that  latitude. 
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VT  had  more  confidence  in  their  measurements  of 
evaporation  than  of  sensible  heat  flux*  Thus,  they  used  the 
Bowen  ratio  to  calculate  sensible  heat  flux  over  the  water, 


ij=^  =  0.66-£-£^  (2.3) 

LHF  1000  e -e.  ^  ’ 


SHF=R  LHF  (2.4) 


where  R  is  the  Bowen  ratio;  SHF  ( cal /cm^ month)  is  the 
sensible  heat  flux;  p  (mb)  is  the  atmospheric  pressure;  Tg 
(C)  and  Tg  (C)  are  the  temperatures  of  the  air  and  sea 
surface . 

The  VT  data  for  wind  speed  were  determined  based  on  ten 
years  of  synoptic  data  at  ship  M  in  the  Norwegian-Barents 
Sea.  Mean  wind  speed  was  based  on  an  evaluation  of  wind 
roses  published  by  the  British  Meteorological  Office  (1959). 
Due  to  the  paucity  of  data,  wind  speed  maps  were  prepared  on 
a  seasonal  basis  and  monthly  values  were  extrapolated.  Air 
and  sea  surface  temperatures  were  obtained  from  the  British 
Meteorological  Office  (1959).  Relative  humidity,  which  is 
needed  to  compute  e^,  was  taken  from  monthly  mean  maps.  The 
relative  humidity  maps  were  produced  based  on  data  from  ship 
M,  small  islands  and  exposed  coastal  stations. 

The  VT  flux  data  are  provided  for  the  Norwegian  and 
Greenland  Seas  in  five  degree  sections  of  latitude,  and  no 
longitudinal  distinctions  are  made.  Data  were  given 
monthly  in  units  of  cal/cm^  month.  Also,  the  total  heat 
fluxes  here  are  computed  by  adding  the  sensible  and  latent 
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fluxes.  All  heat  fluxes  shown  in  this  thesis  are  in  units 
of  W/m^.  To  convert  to  W/m^,  cal/cm^  month  is  multiplied  by 
0.01614.  Tables  2.1,  2.2  and  2.3  are  VT  heat  flux  data  used 
for  this  study. 


Month 

LHF 

SHF 

THF 

January 

59 

142 

201 

February 

69 

171 

240 

March 

89 

206 

295 

Table  2.1  Latent,  sensible  and  total  turbulent  heat  flux 
(W/m^)  over  the  Norwegian-Barents  Sea  from  SON  to  75N  (after 
Vowinckel  and  Taylor,  1964)  . 


Month 

LHF 

SHF 

THF 

January 

73 

45 

118 

February 

64 

60 

124 

March 

90 

86 

176 

Table  2.2  Latent,  sensible  and  total  turbulent  heat  flux 
(W/m^)  over  the  Norwegian-Barents  Sea  from  75N  to  7 ON  (after 
Vowinckel  and  Taylor,  1964)  . 
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Month 

LHF 

SHF 

THF 

January 

113 

44 

157 

February 

84 

47 

131 

March 

95 

54 

149 

November 

67 

32 

99 

Table  2.3  Latent,  sensible  and  total  turbulent  heat  flux 
(W/m^)  over  the  Norwegian-Barents  Sea  from  7 ON  to  65N  (after 
Vowinckel  and  Taylor,  1964) . 


B .  GORSHKOV 

The  World  Ocean  Atlas  edited  by  Gorshkov  (1983), 
contains  annual  and  monthly  maps  of  the  various  flux 
components .  The  atlas  is  in  Russian  with  a  short  summary  in 
English.  The  Russian  headings,  captions  and  other  pertinent 
data  were  translated  by  Capt .  Mary  Lee,  USAF. 

The  main  sources  of  information  for  compiling  the  maps 
were  observations  at  111  shore  stations  from  1936  to  1972. 
The  methodology  for  flux  calculations  was  not  stated. 

The  data  were  provided  as  time  series  depicting  annual 
rates  of  heat  flux  for  twelve  locations  in  kcal/cm^  year. 
Station  10  was  in  the  Norwegian-Barents  Sea,  station  11  was 
in  the  Greenland  Sea  and  station  12  in  the  southern 
Norwegian  Sea  (refer  to  Figure  2.1  for  locations;  all 
figures  are  at  the  end  of  chapters) .  To  convert  to  W/m^, 
kcal/cm^  month  was  multiplied  by  16.14.  Tables  2.4,  2.5  and 
2.6  contain  heat  flux  data  germane  to  this  study. 
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Month 

LHF 

SHF 

THF 

January 

126 

116 

242 

February 

97 

97 

194 

March 

73 

73 

146 

November 

105 

81 

186 

Table  2.4  Latent,  sensible  and  total  turbulent  heat  flux 
(W/m^)  at  station  10  near  72N/020E  over  the  Norwegian- 
Barents  Sea  (after  Gorshkov,  1983) . 


Month 

LHF 

SHF 

THF 

January 

129 

190 

319 

February 

112 

161 

273 

March 

100 

^145 

^245 

Table  2.5  Latent,  sensible  and  total  turbulent  heat  flux 
(W/m  )  at  station  11  near  78N/010E  over  the  Greenland  Sea 
(after  Gorshkov,  1983) . 


Month 

LHF 

SHF 

THF 

January 

129 

186 

315 

February 

113 

161 

274 

March 

97 

65 

162 

November 

89 

61 

150 

Table  2.6  Latent,  sensible  and  total  turbulent  heat  flux 
(W/m  )  at  station  12  near  64N/005E  in  the  Norwegian  Sea 
(after  Gorshkov,  1983). 
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C .  BUNKER 


Bunker  (1976)  computed  energy  flux  over  much  of  the 
world,  including  one  station  in  the  northern  Norwegian  and 
Barents  Seas  that  was  useful  for  this  thesis  (see  Figure 
2.1). 

Bunker  calculated  heat  flux  based  on  the  bulk 
aerodynamic  method.  His  method  used  exchange  coefficients 

vary  with  wind  speed  and  stability.  His  equations  were 
in  the  form 


LHF=  p  ( Os-Qto  ^  ^10  (2.5) 


SHF=  p  ( Ts- )  C/io  (2.6) 


where  p  is  atmospheric  density;  and  Cj,  are  the  variable 
exchange  coefficients  for  water  vapor  and  sensible  heat ;  Cp 
is  the  specific  heat  constant;  Qg  and  Qjq  are  the  averages  of 
the  mixing  ratio  of  air  in  contact  with  the  surface  and  at 
10  m;  Ujo  is  the  average  wind  speed  at  10  m;  and  Tg  and  Tjo 
are  average  temperatures  of  the  sea  surface  and  air  at  10  m. 

Bunkers '  values  for  the  exchange  coefficients  were 
obtained  from  tables  calculated  from  experimental  data.  The 
values  of  the  coefficient  for  water  vapor  have  been  applied 
to  the  computations  of  sensible  heat,  making  the 
coefficients  equal. 

Bunkers'  calculations  were  based  on  a  National  Climatic 
Center  data  set,  which  was  collected  for  the  years  1941- 
1972.  Averages  for  each  month  of  each  year  for  an  entire 
Marsden  Square  were  formed  for  fluxes  and  basic 
meteorological  variables .  Monthly  and  annual  averages  for 
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the  32 -year  period  were  formed  for  subdivisions  of  the 
Marsden  Squares .  A  Norwegian-Barents  Sea  area  was  selected 
at  71  N,  17  E  (refer  to  Figure  2.1  for  location).  This  data 
set  was  produced  based  on  1859  observations  over  the  32-year 
period.  Table  2.7  portrays  the  heat  flux  for  this  station. 


Month 

LHF 

SHF 

THF 

January 

110 

95 

205 

February 

100 

90 

190 

March 

100 

80 

180 

November 

80 

140 

Table  2.7  Latent,  sensible  and  total  turbulent  heat  flux 
(W/m^)  over  the  Norwegian  Sea  at  71N/017E  (after  Bunker, 
1976)  . 


D.  HAKKINEN  AND  CAVALIERI 

Hakkinen  and  Cavalieri  (1989)  estimated  oceanic  surface 
heat  fluxes  in  the  Norwegian,  Greenland  and  Barents  Seas 
using  gridded  Fleet  Numerical  Meteorology  and  Oceanography 
Center  (FNMOC)  surface  analysis  of  winds  and  temperatures 
for  one  year.  They  used  the  bulk  method  (Equations  2.2  and 
2.3)  with  the  heat  (C^)  ,  and  moisture  (Cg)  transfer 
coefficients  fixed  at  a  value  1.5  x  10'^  neglecting  any 
stability  dependence.  Density  was  also  fixed  at  1.3  kg  / 
m^ . 

The  data  for  their  calculations  are  based  on  archived 
Navy  Operational  Global  Atmospheric  Prediction  System 
(NOGAPS)  analyses  for  the  period  December  1,  1978  through 
November  30,  1979.  The  gridded  data  has  a  resolution  of  250 
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km.  The  NOGAPS  model  data  includes  sea  level  pressure,  air 
temperature,  vapor  pressure,  winds  and  SST.  All  parameters 
except  SST  are  given  twice  daily.  SST  is  given  once  daily 
and  interpolated  to  12  hour  intervals  for  flux  computations , 
Monthly  mean  contour  plots  of  total  turbulent  (sensible 
and  latent)  heat  fluxes  were  provided  in  W/m^.  Tables  2.8 
and  2.9  contain  total  turbulent  heat  fluxes  (LHF  +  SHF)  as 
interpreted  from  contoured  maps  for  regions  and  times 
corresponding  with  data  analyzed  in  this  thesis. 


Month 

THF 

January 

640 

February 

560 

March 

420 

Table  2.8  Total  turbulent  heat  flux  (W/m^)  over  the 
Greenland  Sea  77.5N/005E  (after  Hakkinen  and  Cavalieri, 
1989)  . 


Month 

THF 

January 

150 

February 

140 

March 

80 

November 

140 

Table  2.9  Total  turbulent  heat  flux  (W/m^)  over  the 
Norwegian  Sea  68N/005E  (after  Hakkinen  and  Cavalieri,  1989). 
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III.  DATA 


A.  COLLECTION  OP  DATA 


This  study  utilizes  data  collected  by  Professor  Kenneth 
Davidson  and  associates  at  the  Naval  Postgraduate  School 
(NPS)  in  the  Greenland,  Norwegian  and  Barents  Seas  during 
several  recent  ship  cruises  (Figure  3.1).  A  Coastal  Climate 
WeatherPak  meteorological  station  onboard  the  Research 
Vessel  (R\V)  Haakon  Mosby  sampled  air  temperature,  wind 
speed  and  direction  relative  to  the  ship's  heading  every 
second  and  relative  humidity  and  atmospheric  pressure  every 
12  seconds.  The  measurement  platform  was  at  15  m  above  the 
sea  surface.  In  addition,  the  ship's  speed  and  heading  were 
recorded,  enabling  the  calculation  of  true  wind  speed  and 
direction.  Sea  Surface  Temperature  (SST)  was  also  measured 
by  the  ship.  All  observations  were  averaged  over  10  minute 
intervals.  Details  of  the  various  ship  cruises  are 
described  in  later  sections. 

The  temperature  data  were  estimated  to  be  accurate  to 
within  1.0  C  and  the  relative  humidity  data  to  within  5%. 
Wind  speeds  were  estimated  to  be  accurate  to  within  0.3  m/s. 
When  winds  were  from  the  stern,  errors  are  probably  greater 
due  to  the  sheltering  effect  caused  by  the  ship's 
superstructure.  The  measurements  of  atmospheric  pressure 
are  accurate  to  within  2  millibars. 

At  times,  data  were  missed  or  not  recorded;  during 
these  cases  the  trends  were  assumed  to  be  linear  and 
intermediate  data  points  were  interpolated. 

A  common  bias  involving  data  collected  by  ships  is 
storm  avoidance.  That  is,  if  a  ship's  crew  expects  a 
location  to  be  stormy  they  will  typically  avoid  that  area. 
This  is  not  so  with  the  cases  studied  here.  In  general,  the 
locations  and  ship  tracks  were  planned  in  advance  and  were 
not  significantly  deviated  from.  The  only  storm  avoidance 
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bias  occurred  during  the  later  part  of  SIZEX  in  January 
1992,  when  an  intense  cyclone  disabled  the  weather 
collection  unit  with  winds  greater  than  28  m/s.  This 
prevented  about  24  hours  of  data  collection  during  a  high 
wind  period.  This  was  estimated  not  to  significantly  affect 
the  overall  heat  flux  averages. 


EXER 

NORCSEX 

CEAREX 

CEAREX 

NORCSEX 

SIZEX 

SIZEX 

Sea 

Norw 

Barents 

Green 

Norw 

Green 

Barents 

Date 

Mar  88 

Feb  89 

Mar  89 

Nov  91 

Jan  92 

Mar  92 

Wind 

7.0 

10.7 

9.1 

9.0 

10.8 

8.7 

m/s 

(3.5) 

(3.7) 

(2.7) 

(4.0) 

(5.1) 

(3.3) 

Temp 

1.2 

-1.4 

-9.0 

5.8 

-5.6 

0.1 

C 

(2.1) 

(3.9) 

(5.3) 

(2.2) 

(3.2) 

(3.4) 

SST 

5.8 

1.8 

1.4 

8.3 

1.8 

1.8 

C 

(1.1) 

(0) 

(0) 

(0.5) 

(2.0) 

(3.0) 

RH 

68.6 

72.0 

72.4 

69.3 

71.7 

81.5 

% 

(11.3) 

(12.9) 

(9.6) 

(13.5) 

(11.9) 

(9.7) 

Press 

1002 

997 

991 

999 

996 

1001 

mb 

(8) 

(8) 

(11) 

(14) 

(12) 

(9) 

Table  3.1  Mean  and  standard  deviations  of  input  data  for 
indicated  exercises. 


B.  NORCSEX  88 

The  Norwegian  Continental  Shelf  Exercise  (NORCSEX)  88 
was  a  European  Research  Satellite  (ERS-1)  pre-launch 
calibration  and  validation  experiment.  The  purpose  of  the 
ship  cruise  was  to  validate  data  measured  by  aircraft  flown 
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with  ERS-1  sensors.  NORCSEX  88  was  conducted  in  the 
Norwegian  Sea  and  centered  around  the  Norwegian  Continental 
Shelf  at  64°30'N  and  009°E  during  March  of  1988.  From  07  to 
18  March,  the  R/V  Haakon  Mosby  collected  meteorological  data 
while  transiting  along  the  west  coast  of  Norway  (see  Figure 
3.1  for  track).  The  time  series  of  these  data  are  shown  in 
Figure  3.2. 

NORCSEX  88  displayed  the  lowest  wind  speeds  of  all  the 
cases;  during  the  exercise  moderate  winds  prevailed  with  a 
mean  of  only  7.0  m/s.  However,  gale  to  storm  force 
conditions  occurred  on  11  March  (Julian  day  71);  as  the  wind 
veered  from  the  south  to  the  north,  the  wind  speed  quickly 
increased  to  more  than  22  m/ s .  Along  with  the  wind  change 
came  a  decrease  in  air  temperature  and  an  increase  in 
relative  humidity.  Throughout  the  measurement  period,  the 
air  temperature  and  sea  surface  temperature  were  fairly 
constant.  SST  had  a  range  of  2  C  to  8  C  and  a  standard 
deviation  of  1.1.  Air  temperature  had  a  range  of  -5  C  to  5 
C  and  the  smallest  standard  deviation  of  all  the  exercises 
at  2.1  C. 

C .  CEAREX 

The  Coordinated  Eastern  Arctic  Exercise  (CEAREX)  was  a 
large  international  interdisciplinary  project  involving 
several  ships,  aircraft  and  ice  camps.  CEAREX  was  conducted 
in  the  vicinity  of  the  Svalbard  Islands  from  September  1988 
to  May  1989.  For  this  study,  continuous  surface 
measurements  collected  by  the  R/V  Haakon  Mosby  from  25  Feb 
to  23  Mar  1989  were  used.  From  25  February  to  09  March  1989 
the  ship  was  in  the  northern  Norwegian  and  Barents  Seas 
(Figure  3.1).  Since  these  tracks  correspond  more  closely  to 
land  and  climatological  stations  in  the  Barents  Sea,  this 
region  will  be  called  only  the  Barents  Sea  to  prevent 
confusion  with  locations  in  the  southern  Norwegian  Sea. 
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After  10  March,  the  ship  transited  on  a  westerly  course  to  a 
location  west  of  Svalbard  and  conducted  north-south  tracks 
along  the  ice  edge.  Since  the  measurements  in  the  Barents 
Sea  differed  markedly  from  those  in  the  Greenland  Sea,  the 
measurements  have  been  divided  into  two  separate  data  sets. 
The  time  series  of  data  in  the  Barents  Sea  is  shown  in 
Figure  3.3,  and  the  time  series  for  data  in  the  Greenland 
Sea  is  shown  in  Figure  3.4. 

There  were  several  cyclone  passages  while  the  ship  was 
located  in  the  Barents  Seas .  The  temperature  tended  to  be 
near  freezing  for  relatively  long  periods  of  time  with  brief 
surges  of  extremely  cold  air  after  cyclone  passages .  The 
average  wind  was  strong  at  10.7  m/ s  and  primarily  from  the 
north.  Sea  surface  temperatures  were  not  available  for 
CEAREX.  Mean  SSTs  of  1.8  C  and  1.4  C  were  chosen  for  CEAREX 
data  in  the  Barents  and  Greenland  Seas,  respectively.  These 
values  were  based  upon  averages  from  SIZEX  data  along 
similar  tracks  and  climotological  data. 

In  the  Greenland  Sea,  the  air  temperatures  were  much 
colder  than  in  the  Norwegian  Sea,  with  mean  temperatures  of 
-9  C  compared  with  -1  C.  After  passage  of  a  cyclone  on 
about  11  March  (Julian  day  69)  the  pressure  field  was 
relatively  flat  with  colder  temperatures  and  lower  relative 
humidity. 

D.  NORCSEX  91 


The  Norwegian  Continental  Shelf  Exercise  (NORCSEX)  91 
was  a  follow-on  to  NORCSEX  88.  Its  purpose  was  a  post 
launch  validation  of  ERS-1  sensors.  NORCSEX  91  was 
conducted  in  the  Norwegian  Sea  from  08  Nov  to  29  Nov  1991. 
During  this  period  meteorological  variables  were  collected 
onboard  the  R/V  Haakon  Mosby.  The  time  series  of  data  are 
shown  in  Figure  3.5. 

NORCSEX  91  data  contained  the  warmest  air  and  sea 
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surface  temperatures  of  all  the  exercises  with  means  of  5.8 
C  and  8.3  C,  respectively.  The  temperatures  were  also  less 
variable  than  those  of  other  exercises  with  standard 
deviations  of  2.2  C  and  0.5  C.  During  the  period  there  were 
several  strong  cyclones  that  passed  over  the  ship.  These 
cyclones  exhibited  higher  winds  but  showed  little 
temperature  change. 

E .  SIZEX 

The  Seasonal  Ice  Zone  Experiment  (SIZEX)  was  conducted 
as  a  part  of  an  ERS-1  ice  Synthetic  Aperture  RADAR  (SAR) 
validation.  SIZEX  was  held  in  the  Greenland  Sea  during 
1992.  The  experiment  has  been  divided  into  two  periods,  one 
in  January  and  one  in  March.  Time  series  of  data  are  shown 
in  figures  3.6  and  3.7. 

In  January,  the  data  were  influenced  by  three  prominent 
cyclone  passages.  The  cyclones  were  accompanied  by  strong 
winds  from  the  south  to  southeast.  The  final  cyclone, 
occurring  on  about  13  January  (Julian  day  13),  was  storm 
force  with  wind  speeds  approaching  30  m/s.  The  storm  was 
strong  enough  to  knocJc  out  the  weather  collection  station. 

In  February  and  March,  there  were  no  intense  cyclones 
but  periods  with  significant  winds  between  15-20  m/s. 
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julian  date 


igure  3.3  Time  series  of  wind  speed,  wind  direction,  sea 
urface  temperature  (dotted)  and  air  temperature  (solid), 
elative  humidity  and  pressure  during  CSAREX  in  the  Barents 
ea,  February -March  1989. 
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Figure  3.4  Time  series  of  wind  speed,  wind  direction,  se 
surface  temperature  (dotted)  and  air  temperature  (solid) 
relative  hiamidity  and  pressure  during  CEAREX  in  the 
Greenland  Sea,  March  1989. 
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Figure  3.5  Time  series  of  wind  speed,  wind  direction,  s 
surface  temperature  (dotted)  and  air  temperature  (solid) 
relative  humidity  and  pressure  during  N0RCSEX91  in  the 
Norwegian  Sea,  November  1991. 
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Figure  3.6  Time  series  of  wind  speed,  wind  direction,  sea 
surface  temperature  (dotted)  and  air  temperature  (solid) , 
relative  humidity  and  pressure  during  SIZEX  in  the  Greenland 
Sea,  January  1992. 
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Figure  3.7  Time  series  of  wind  speed,  wind  direction,  sea 
surface  temperature  (dotted)  and  air  temperature  (solid) , 
relative  humidity  and  pressure  during  SIZEX  in  the  Barents 
Sea,  March  1992. 


IV.  HEAT  FLUX  COMPUTATIONS  USING  THE  BULK  METHOD 


A.  THEORY 


The  primary  goal  of  this  thesis  is  to  quantify  the 
turbulent  heat  flux  between  the  ocean  and  atmosphere  in  the 
Greenland,  Norwegian  and  Barents  Sea  regions.  This  goal  is 
accomplished  through  use  of  the  bulk  aerodynamic  method 
which  is  based  on  Monin-Obukhov  (M-0)  similarity  theory. 

The  bulk  method  has  proven  to  be  an  accurate  method  for 
computing  heat  fluxes  between  the  ocean  and  atmosphere 
(Blanc,  1985) .  This  method  is  the  most  practical  for 
shipboard-based  measurements,  requiring  only  routine 
meteorological  measurements  such  as  wind  speed,  air 
temperature,  sea  surface  temperature,  relative  humidity  and 
atmospheric  pressure.  The  bulk  method  allows  calculations 
to  be  made  based  on  measurements  made  from  ships  at  sea, 
which  provide  a  substantial  portion  of  the  data  available 
over  the  oceans . 

Monin-Obukhov  similarity  theory  assumes  the  gradients 
of  the  meteorological  variables  depend  on  surface  fluxes, 
height  and  stability  only.  For  instance,  the  wind  is 
modeled  using  the  log  wind  profile  modified  for  stability 
effects.  Therefore,  if  the  wind  speed  is  known  at  any 
height  within  the  surface  layer,  one  can  estimate  the  wind 
speed  at  any  other  height  within  the  surface  layer  (Equation 
4.1) .  Likewise,  the  profiles  of  temperature  and  moisture 
can  be  used  to  estimate  sensible  and  latent  heat  fluxes . 

The  profiles  for  temperature  and  humidity  are  known  and  can 
be  found  by  Equation  4.2, 
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(4. 1) 


du 

dz 


kz 


4) 


Af 


dr 

dz 


(4.2) 


Where  u(in/s),  T  (C)  anci  z  (m)  are  wind  speed,  potential 
temperature  and  observation  height;  u.  (m/s)  is  a  velocity 
scale  called  the  friction  velocity,  T.  (C)  is  a  surface 
layer  temperature  scale;  k  is  Von  Karman's  constant;  and 
and  dimensionless  wind  shear  and  temperature 

functions  of  stability  (z/L) . 

B .  BULK  FORMULA 


The  turbulent  heat  fluxes  at  the  surface  can  be 
expressed  in  the  form 

SHF  =  pCp  w'T'=  -pCpU,  r.  (4.3) 


LHF=pL^w'g'=-pL^u,g,  (4.4) 

where  SHF  (W/m^  )  and  LHF  (W/m^  )  are  sensible  and  latent 
heat  flux,  with  the  sign  convention  that  a  positive  value 
indicates  a  heat  flux  from  the  ocean  to  the  atmosphere  and  a 
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negative  indicates  heat  flux  from  the  air  to  the  ocean;  The 
variable  p  (kg/m^)is  atmospheric  density,  based  on  surface 
temperature  and  pressure  and  calculated  using  the  ideal  gas 
law;  Cp  (1004  J/kg  K)  is  the  specific  heat  of  dry  air; 

(250000  J/kg) is  the  latent  heat  of  vaporization  at  0°C. 

The  parameter  T.  (C)  is  a  temperature  scale  and  q,  (g 
water/g  air)  is  a  humidity  scale.  The  scaling  paramters  are 
calculated  using  the  following  equations 


u 


* 


uk 

)-T, 

^0 


(4.5) 


(4.6) 


C[,= 


ig^o-gjk 


(4.7) 


where  Tjo  (C)  and  Tg  (C)  are  potential  temperature  at  10  m 
and  at  the  surface;  q^  (g/kg)and  q^  (g/kg)  are  the  mixing 
ratios  at  10  m  and  at  the  surface;  Zq  is  roughness  length, 
which  is  the  sum  of  Charnock's  roughness  length  and  the 
roughness  length  for  a  smooth  surface  and  is  defined  by 
Equation  4.8;  Zq^  and  z  oq  are  the  thermal  and  moisture 
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roughness  lengths  based  on  a  constant  heat  exchange 
coefficient  (C^n  and  )  and  defined  by  Equations  4.9  and 
4.10; 

Zo  =  (0.011u.V9.8) +(  (.11)  (1.4x10-5) /u,)  (4.8) 


=  exp 


Ctn  ln(z/2o) 


(4.9) 


.„,=  10exp— 4.10 


and  are  the  integrated  stability  functions,  for 
unstable  stratification  and  defined  by  Equations  4.11  and 
4.12  . 

In  [  (1  +x)  /2]  +ln  [  (1  +x^)  /2]  -2  arctan(x)  +it/2  (4.11) 

7t  =  2  In  [  (l+x2) /2]  (4.12) 


x=(l-i6z/L) -25  (4.13) 
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and  for  stable  stratification 


Yj„  =  -5  z/L  (4.14) 
Tt  =  -5  z/L  (4.15) 


L  is  the  Monin-Obukhov  stability  length  defined  by 
Equation  4.16,  and  z/L  is  the  stability;  T^,  is  the  buoyancy 
flux  scaling  parameter  determined  by  Equation  4.17, 


L  = 


(4.16) 


T  -  (r,-r,,)k 

ln(^/zo)-Yj. 


(4.17) 


where  (C)  and  (C)  are  the  virtual  potential 
temperature  at  the  measurement  height  and  the  surface. 

The  constants  in  the  Charnock  relationship  (Equation 
4.8)  and  the  integrated  stability  functions  are  from  Smith 
(1988)  . 

The  advantage  of  using  the  bulk  method  is  that  the  only 
P^^^nieters  required  are  the  mean  wind  speed,  the  mean  air 
temperature  and  the  mean  humidity  at  one  level  each  in  the 
surface  layer,  and  the  mean  surface  temperature.  High 

multi-level  measurements  are  not  required  and 
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the  bulk  method  is  less  affected  by  ship  effects  than  are 
other  flux  calculation  methods. 

For  the  unstable  cases  (cold  air  over  warm  water)  and 
stable  cases  (warm  air  over  cold  water)  'Ft  is  a  function  of 
the  stability  parameter  (z/L) .  The  surface  layer  over  open 
water  in  the  Arctic  is  typically  unstable.  Therefore,  z/L 
is  negative,  and  the  sensible  and  latent  heat  fluxes  are 
enhanced  by  buoyancy  effects.  However,  if  the  surface  layer 
is  stable,  which  is  rare  but  does  occur,  then  z/L  is 
positive.  Thus,  'Ft  is  negative  and  u.,  T*  and  q.  would  all 
be  smaller  in  magnitude;  this  implies  the  sensible  and 
latent  heat  fluxes  would  also  be  smaller  in  magnitude. 

C.  SOLVING  FOR  HEAT  FLUX  BY  ITERATION 

Determining  u,,  T*,  and  g.  requires  knowledge  of  L, 
which  in  turn  is  a  function  of  u.,  T*  and  q, .  Therefore,  an 
iterative  solution  method  is  needed. 

The  first  step  in  the  iterative  process  is  to  make  an 
initial  guess  at  L,  Zq,  and  u, .  Then  Equations  4.5,  4.6, 

4.7,  4.8  and  4.16  are  repeatedly  solved  for  u,,  Zq  and  L 
until  the  solution  for  u*  converges.  However,  if  the  wind 
speed  is  zero  or  the  atmosphere  is  too  stable,  the  solution 
will  not  converge.  In  this  case  the  atmosphere  is  non- 
turbulent  and  u,,  T.,  q.  and  z/L  are  set  equal  to  zero.  The 
Matlab  program  used  in  calculating  heat  fluxes  can  be  found 
in  the  Appendix . 

D.  LIMITATIONS  OF  THE  BULK  METHOD 

The  scheme  shown  above  is  one  of  many  bulk  transfer 
coefficient  schemes .  Each  of  the  different  schemes  contains 
stent  coefficients  of  drag,  heat  and  humidity  and  varies 
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in  degree  of  complexity.  Some  include  stability  dependence 
while  others  ignore  stability  or  assume  a  constant  stability 
effect.  Additionally,  each  method  has  a  different  range  of 
wind  speeds  and  air  sea  temperature  differences  under  which 
the  method  is  considered  valid.  Smith  (1988)  listed  several 
limitations  to  his  technique,  which  are  summarized  below. 

1.  High  Wind  Speed 

There  are  few  reliable  measurements  of  wind  stress  and 
heat  flux  above  26  m/s.  Thus,  application  of  the  method  to 
wind  speeds  above  26  m/ s  would  be  an  extrapolation  of  the 
lower  wind  speed  data.  Fortunately,  winds  above  the  26  m/s 
threshold  are  absent  in  all  but  the  SIZEX  data  set  and  are 
rare  there. 

2.  Extremely  Stable  Stratification 

If  the  atmosphere  is  too  stable  then  the  iterative 
method  will  fail  to  produce  a  convergent  solution.  Here  the 
program  will  give  a  value  of  zero  for  u*,  T*,  q*  and  z/L. 

This  implies  the  fluxes  are  zero  and  actual  values  will  be 
very  close  to  zero.  These  highly  stable  conditions  would 
occur  if  the  air  temperature  was  very  warm  and  the  sea 
surface  very  cold,  but  such  occurrences  are  rare. 

3.  Extremely  Unstable  Stratification 

If  the  wind  speed  were  close  to  zero  then  the  iterative 
method  would  fail  to  produce  a  convergent  solution.  In  this 
case  the  program  will  give  a  value  of  zero  for  u*,  which 
implies  the  fluxes  are  zero.  In  fact,  there  may  be  fluxes 
due  to  convection.  Free  convective  scaling  was  not  included 
in  the  bulk  formula,  but  this  was  judged  to  not  degrade  the 
results  since  calm  conditions  and  non- convergent  solutions 
are  very  rare  or  nonexistent  in  the  data. 
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V.  RESULTS  AND  COMPARISONS 


A.  RESULTS  OF  HEAT  FLUX  COMPUTATIONS 


Turbulent  heat  fluxes  were  computed  based  on  the  data 
presented  in  Chapter  III  and  the  method  discussed  in  Chapter 
IV.  Overall,  the  average  measured  turbulent  heat  fluxes 
(THF)  in  this  thesis  (Table  5.1)  compared  favorably  with 
those  given  by  the  climatological  studies  mentioned 
previously . 


Exercise 

Date 

Region 

LHF 

SHF 

THF 

CEAREX 

Mar  89 

Greenland  Sea 

109 

140 

249 

SIZEX 

Jan  92 

Greenland  Sea 

111 

105 

216 

NORCSEX 

Mar  88 

Norwegian  Sea 

87 

44 

131 

NORCSEX 

Nov  91 

Norwegian  Sea 

96 

30 

126 

CEAREX 

Feb  89 

Barents  Sea 

73 

45 

118 

SIZEX 

Mar  92 

Barents  Sea 

44 

20 

64 

Table  5.1  Mean  value  of  latent,  sensible  and  total 
turbulent  heat  flux  (W/m^)  . 
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Heat  flux  computations  for  all  the  NFS  ship  cruises 
occurred  during  the  Arctic  winter  (November  -  March) . 
Analysis  of  the  climatological  data  of  Bunker  (1976)  and 
Gorshkov  (1988),  showed  peak  heat  loss  (positive  or  upward 
heat  flux)  for  this  region  in  December  and  January  with  only 
slightly  less  heat  loss  in  the  preceding  and  following 
months.  Therefore,  the  time  difference  in  the  observations 
does  not  significantly  hinder  comparisons  between  exercises 
that  occurred  in  different  winter  months,  and  the 
differences  in  heat  flux  values  displayed  in  Table  5.1  are 
primarily  due  to  the  differenecs  in  the  ship's  location  or 
differences  between  years. 

According  to  the  data,  the  largest  mean  THF,  near  250 
W/m^  was  seen  in  the  Greenland  Sea.  The  Norwegian  Sea  had 
a  mean  THF  of  130  W/m^,  whereas  the  Barents  Sea  had  the 
lowest  mean  THF. 

B.  COMPARISON  WITH  PREVIOUS  STUDIES 

In  this  section,  the  calculations  of  heat  flux  based  on 
the  NFS  data  are  compared  with  previously  published  values. 
Accepted  faults  with  comparisons  include:  (1)  Slightly 
differing  spatial  regimes:  the  mean  value  of  the  total 
turbulent  heat  flux  measurement  is  taken  along  a  ship  track, 
whereas  in  the  climatological  studies  locations  are  fixed  or 
averaged  over  large  areas  and  HC  used  grid  point  data.  (2) 
Different  time  scales:  record  length  was  two  to  three  weeks, 
whereas  the  climatological  studies  have  tens  of  years  of 
data,  and  HC  data  comes  entirely  from  one  year  (1979)  . 
However,  the  bulk  method  used  here,  is  similar  to  Bunkers' 
method,  and  is  considered  more  advanced  than  some  earlier 
methods.  For  example,  the  method  used  in  this  thesis 
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includes  flux  computations  that  account  for  variable 
stability  (z/L),  roughness  (Zq)  and  atmospheric  density  (p) . 

In  the  following  subsections,  the  results  of  each 
exercise's  heat  flux  computations  are  reviewed  and  compared, 
by  region,  to  the  results  of  the  earlier  studies. 

1 .  Greenland  Sea 

During  SIZEX,  the  mean  LHF  was  111  W/m^  and  the  mean 
SHF  was  105  W/m^.  This  exercise  had  some  of  the  largest  10 
minute  averaged  heat  fluxes  of  all  the  cases  (Figure  5.1). 
These  large  heat  fluxes  are  a  result  of  the  stronger  winds 
and  strong  to  moderate  air  sea  temperature  differences 
(ASTD)  encountered  in  the  Greenland  Sea  during  SIZEX. 
Histograms  of  THF  for  all  exercises  were  plotted  with  50 
equally  spaced  bins  between  the  maximum  and  minimum  values 
of  THF  for  that  exercise.  The  histogram  for  this  case 
(Figure  5.2)  showed  a  maximum  occurrence  near  the  median  of 
210  W/m^  and  a  long  tail  to  the  right  (above  400  W/m^’ 
corresponding  to  a  high  wind  event  on  Julian  day  14.  A 
false  local  maximum  occurring  at  100  W/m^  is  the  result  of  a 
data  gap  of  about  six  hours  on  Julian  day  13 . 

During  CEAREX,  the  heat  fluxes  (Figure  5.3)  were  also 
relatively  large  compared  with  other  NFS  exercises .  The 
mean  LHF  was  109  W/m^  and  the  mean  SHF  was  140  W/m^.  As  in 
SIZEX,  the  large  heat  fluxes  can  be  attributed  to  large 
ASTD,  large  moisture  gradients  and  moderate  winds.  The 
large  ASTD  values  are  a  result  of  the  cold  air  advection  off 
the  ice  (northerly  winds).  Mean  winds  (9.1  m/s)  for  the 
Greenland  Sea  region  during  CEAREX  are  similar  to  other 
exercises  and  had  a  small  standard  deviation  (2.7  m/s) 
compared  with  other  exercises.  Nevertheless,  the  heat 
fluxes  are  larger  in  the  Greenland  Sea  primarily  due  to  the 
colder  air  over  the  Greenland  Sea.  The  histogram  of  the 
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heat  flux  (Figure  5.4)  has  at  least  three  peaks.  The 
different  peaks  are  the  result  of  changes  in  synoptic 
conditions.  The  largest  number  of  occurrences  is  at  330 
VI/w?  and  is  related  to  the  cold  air  outbreak  on  Julian  days 
71  through  75.  During  this  period  the  highest  heat  fluxes 
were  caused  by  cold  air  advection,  and  occurred  after  the 
winds  backed  to  the  north  and  result  in  a  highly  unstable 
surface  layer  (very  cold  air  over  warmer  water) .  The  second 
peak  is  near  200  W/m^  and  is  the  result  of  smaller  winds  and 
ASTD  after  Julian  day  75"  The  third  local  maximum  at  60 
W/m^  occurred  during  a  pressure  minimum,  and  was  the  result 
of  the  southerly  winds  that  gave  an  unstable  surface  layer 
(cold  air  over  warmer  water) . 

The  VT  heat  flux  data  compared  well  with  the  January 
SIZEX  data,  but  were  20%  larger  than  the  March  CEAREX 
observations  of  heat  flux  (Tables  5.2  and  5.3).  The 
Gorshkov  data  showed  mixed  results.  Their  values  for  heat 
flux  during  SIZEX  were  within  10%  and  during  CEAREX  were  70% 
less  than  computed  for  the  NFS  data  set.  The  HC  heat  fluxes 
were  two  to  three  times  larger  than  any  other  computed  heat 
flux  during  January,  and  were  20%  larger  than  the  NFS  values 
during  CEAREX,  which  was  in  March. 
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Exercise 

LHF 

SHF 

THF 

SIZEX  GS 

111 

105 

216 

VT 

59 

142 

201 

Gorshkov 

126 

116 

242 

HC 

NA 

NA 

640 

Table  5.2  Latent,  sensible  and  total  turbulent  heat  flux 
(W/m^)  over  the  Greenland  Sea  in  January. 


Exercise 

LHF 

SHF 

THF 

CEAREX  GS 

109 

140 

249 

VT 

89 

206 

295 

Gorshkov 

73 

73 

146 

HC 

NA 

NA 

300 

Table 

(W/m^) 


d.3  Latent,  sensible  and  total  turbulent  heat  flux 
over  the  Greenland  Sea  in  March. 


2 .  Norwegian  Sea 

NORCSEX  88  was  conducted  in  March  1988,  off  the 
southwest  coast  of  Norway .  Heat  fluxes  for  this  exercise 
were  similar  to  previous  results  for  the  area  (Table  5.4). 
The  mean  LHF  was  87  W/m^  and  the  mean  SHF  was  44  W/m^. 

The  heat  flux  (Figure  5.5)  exhibits  a  synoptic 
influence,  which  should  be  expected  in  a  storiry  region  of 
the  Norwegian  Sea.  For  example,  before  the  passage  of  a 
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series  of  surface  cyclones  on  Julian  day  71  the  THF  is  small 
(  50  W/m^  )  due  to  low  winds  and  a  small  ASTD.  Afterwards, 
the  wind  speed  and  ASTD  increased  causing  the  THF  to  climb 
above  2  00  W/m^  and  then  decreased  to  around  the  mean  of  13  0 
W/m^  after  the  winds  died  down.  Also,  the  SHF  time  series 
is  strikingly  similar  to  the  time  series  of  wind  speed  (as 
seen  in  Chapter  III,  Figure  3.1) .  The  histogram  of  heat 
flux  (Figure  5.6)  is  nearly  symmetric  with  a  sharp  peak  near 
the  median  of  120  W/m^.  The  tail  of  fluxes  above  250  W/m^ 
represents  a  mesoscale  gale  force  storm  observed  for  several 
hours  on  Julian  day  71. 

NORCSEX  91  was  conducted  in  the  southern  Norwegian  Sea 
in  November  1991.  The  THF  values  for  this  exercise  (Figure 
5.7)  were  similiar  to  those  of  NORCSEX  88.  During  NORCSEX 
91,  the  SHF  was  low  with  a  mean  30  W/m^.  The  LHF  was 
larger,  with  a  mean  value  of  96  W/m^.  NORCSEX  91 
experienced  the  passage  of  several  cyclones  with  associated 
wind  variations.  The  heat  fluxes  modulated  with  these  wind 
changes.  Also,  apparent  is  the  effect  of  the  ASTD.  During 
the  earlier  periods  with  larger  ASTD  the  fluxes  are  larger 
but  after  Julian  day  328  when  the  ASTD  became  smaller  the 
flux  decreased.  The  histogram  of  THF  (Figure  5.8)  has  three 
modes.  The  first  mode,  near  0  W/m^  was  a  result  of  the 
small  or  zero  ASTD  that  occurred  from  Julian  day  333  to  335. 
Modes  two  and  three  occurring  at  90  W/m^  and  170  W/m^  are 
the  result  of  a  bimodal  distribution  about  the  median  of  125 
W/m^- 

The  differing  data  sets  compare  well  within  this  region 
as  all  the  data  sets  are  within  20%  -  30%  of  the  mean 
observed  THF  (Tables  5.4  and  5.5). 
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Exercise 

LHF 

SHF 

THF 

NORCSEX  88 

87 

44 

131 

VT 

95 

54 

149 

Gorshkov 

97 

65 

162 

HC 

NA 

NA 

100 

Table  5.4  Latent,  sensible  and  total  turbulent  heat  flux 
(W/m^)  over  the  Norwegian  Sea  in  March. 


Exercise 

LHF 

SHF 

THF 

NORCSEX  91 

96 

30 

126 

VT 

67 

32 

99 

Gorshkov 

89 

61 

150 

HC 

NA 

NA 

160 

Table  5.5  Latent,  sensible  and  total  turbulent  heat  flux 
(W/m^)  over  the  Norwegian  Sea  in  November. 


3 .  Barents  Sea 

A  segment  of  CEAREX  was  conducted  in  the  Barents  Sea  in 
February  and  early-March  1989.  During  CEAREX  in  the  Barents 
Sea  the  mean  LHF  was  73  W/m^  and  the  mean  SHF  was  45  W/m^. 
Several  cyclones  seen  in  the  data  are  apparent  in  the  time 
series  of  heat  flux  as  peaks  at  intervals  of  two  to  three 
days  (Figure  5.9).  The  histogram  of  THF  (Figure  5.10)  is 
skewed  to  the  right  and  has  some  high  peaks  at  80  W/m^ . 
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Also  noticeable  are  the  negative  fluxes  (heat  flux  into  the 
ocean)  resulting  from  the  surge  of  warm  southerly  air  and 
stable  conditions  around  Julian  day  66. 

SIZEX  was  conducted  in  the  Barents  Sea  in  March.  The 
THF  values  for  March  in  the  Barents  Sea  were  small,  near  60 
W/m^.  During  SIZEX,  the  LHF  was  44  W/m^  and  the  SHF  was  20 
W/m^.  The  heat  fluxes  are  smaller  here  due  to  smaller  ASTD 
(Figure  5.11).  The  histogram  of  THF  (Figure  5.12)  is 
symmetric  with  a  median  of  near  60  W/m^. 

The  heat  fluxes  of  VT  and  Bunker  are  up  to  two  times 
greater  than  observed  by  NFS  ship  data  in  the  Barents  Sea 
(Tables  5.6  and  5.7).  However,  the  Gorshkov  and  HC  heat 
fluxes  are  three  to  four  times  larger  than  the  measured 
values.  The  values  of  heat  flux  in  the  NFS  data  are 
strongly  influenced  by  the  lack  of  SST  data  and  the 
assumption  of  a  constant  1.8  C  SST  for  CEAREX.  For  example, 
if  the  average  values  for  all  other  directly  measured 
quantities  (wind  speed,  air  temp,  rh,  p)  were  unchanged  and 
the  average  SST  was  assumed  to  be  0  C  vice  1.8  C,  the  THF 
would  become  60  W/m^,  nearly  half.  If  the  average  SST  were 
5  C,  the  THF  would  become  230  W/m^  or  nearly  double  the 
estimated  value.  Normally  SST  would  not  be  so  important. 
However  in  this  instance  it  is  critical  because  the  assumed 
SST  (1.8  C)  and  the  actual  mean  air  temps  (-1.4  C)  are 
similar . 
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Exercise 

LHF 

SHF 

THF 

CEAREXBS 

73 

45 

118 

VT 

64 

60 

124 

Gorshkov 

112 

161 

273 

Bunker 

100 

90 

190 

HC 

NA 

NA 

400 

Table  5.6  Latent,  sensible  and  total  turbulent  heat  flux 
(W/m^)  over  the  Barents  Sea  in  February. 


Exercise 

LHF 

SHF 

THF 

SIZEX  BS 

44 

20 

64 

VT 

90 

86 

176 

Gorshkov 

100 

145 

245 

Bunker 

100 

80 

180 

HC 

NA 

NA 

240 

Table  5.7  Latent,  sensible  and  total  turbulent  heat  flux 
(W/m^)  over  the  Barents  Sea  in  March. 
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C.  LIMITATIONS  OF  PREVIOUS  STUDIES 


1.  Time  Averaging  Methods 

Hanawa  and  Toba  (1987)  (hereafter  HT)  made  comparisons 
of  data  averaged  over  various  time  periods  from  one  day  to  a 
month  using  the  sampling  method  (SM)  and  the  scalar 
averaging  method  (SAM) .  The  sampling  method  (SM)  is  a 
simple  average  of  fluxes  computed  from  hourly  (or  every 
observation)  values  of  the  bulk  meteorological  parameters 
(wind,  temperature  and  humidity) .  The  scalar  averaging 
method  computes  the  average  flux  from  average  values  of  the 
bulk  data,  thus  ignoring  correlations  between  variables. 

Both  methods  were  compared  with  data  taken  every  three  hours 
from  a  mid-latitude  Oceanographic  Weather  Station  (OWS-T) 

(29N  135E)  from  June  1950  to  November  1953.  HT  showed  that 
the  SHF  was  almost  the  same  regardless  of  method  and 
averaging  time.  Also,  LHF  values  computed  by  the  SAM  are 
about  105%  of  those  by  the  SM  for  averaging  times  from  three 
days  to  one  month. 

Ledvina  et  al .  (1993)  made  comparisons  similar  to  HT, 

using  averaging  times  from  2-72  hours.  Their  ratios  vary 
little  for  times  exceeding  36  hours,  thereby  eliminating  the 
need  to  extend  the  averaging  period  beyond  72  hours.  Their 
data  was  in  the  equatorial  Pacific  near  OOON  145E  from  17 
February  to  10  March  1990.  Like  HT,  they  showed  the  SAM 
overestimated  LHF  with  increasing  averaging  period.  Yet 
with  results  drastically  different  from  HT,  they  showed  the 
SAM  to  underestimate  SHF  by  40%.  Much  of  this  difference 
may  be  related  to  convection  in  the  tropics,  thereby 
demonstrating  the  importance  of  regional  influences  on  heat 
fluxes . 

The  purpose  of  this  subsection  is  to  extend  the  work 
completed  by  HT  and  Ledvina  et  al .  (1993)  to  the  polar 
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latitudes  with  recent  data,  thus  providing  a  basis  for 
determining  the  accuracy  of  heat  fluxes  computed  via  various 
averaging  techniques  in  the  polar  seas.  Ideally,  time 
averages  are  desired  that  are  long  enough  to  include  all  low 
frequency  effects  (synoptic)  and  sampled  fast  enough  to 

capture  all  high  frequency  effects  (Kaimal  and  Finnigan, 

1994)  . 

Esbensen  and  Reynolds  (1981)  define  the  sampling  method 
for  calculating  SHF  and  LHF  in  the  following  equations 

WiF  =  pCpC^C/AT  (5.1) 

LHF  =  pLyC^lFKg  (5.2) 


where  averaging  is  indicated  by  the  overbar.  The  key 
feature  of  this  method  is  that  heat  fluxes  are  calculated 
for  every  measurement.  The  measurement  interval  was  10 
minutes  in  the  NFS  data  set,  generally  3  hours  in  Bunker's 
data  to  12  hours  in  the  HC  model .  The  fluxes  are  then 
summed  and  then  averaged. 

The  scalar  averaging  method  (classical  method)  for 
computing  heat  flux  defines  the  flux  in  the  following  forms 

SHF  =  pCpC^U'Kf  (5.3) 


LHF=pLyC^U'^  (5.4) 
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where  the  average  transfer  coefficients  of  heat  and  humidity 
are  computed  empirically  from  the  sample  means  of 
temperature  difference,  humidity  difference  and  wind  speed. 
The  difference  of  the  classical  method  (SAM)  from  the 
sampling  method  is  the  measurements  are  averaged,  usually 
over  the  entire  month  and  then  the  heat  fluxes  are  computed 
based  on  the  averages  of  the  measurements . 

Lumley  and  Panofsky  (1964)  suggest  that  a  reasonable 
averaging  time  for  average  mean  horizontal  winds  (and 
temperature)  would  be  30  minutes,  corresponding  to  the 
passage  of  two  to  three  large  convective  cells  through  the 
depth  of  the  convective  boundary  layer.  However,  if 
averaging  periods  are  increased  to  about  one  hour, 
nonstationarity  in  the  form  of  diurnal  variations  will 
occur.  With  averaging  periods  of  one  hour  or  greater, 
errors  of  5%  and  10-50%  for  T.  and  u.  result. 

Bunker  (1976)  and  HC  used  the  sampling  method  when 
calculating  their  heat  fluxes;  VT  and  Gorshkov  (1988)  used 
the  classical  approach  or  scalar  averaging  method.  Finally, 
data  from  this  study  are  used  to  compare  the  sampling  method 
with  the  classical  method  in  the  Greenland,  Norwegian  and 
Barents  Seas. 

The  approach  of  Bunker  (1976)  was  to  include  every  ship 
or  buoy  report  and  compute  heat  fluxes  and  then  average  the 
fluxes  over  the  month.  This  method  was  better  than  the 
classical  approach  in  that  it  would  not  result  in  errors  due 
to  diurnal  fluctuations.  The  classical  approach  used  by 
Gorshkov  (1988)  and  VT  also  has  it  advantages.  Since  the 
covariance  and  correlation  of  UAT  and  UAq  are  small  at 
shorter  time  intervals  and  increase  for  longer  periods,  the 
classical  method  may  provide  a  better  heat  flux  estimate 
when  the  data  set  is  small  (Ebsensen  and  Reynolds,  1981) . 

The  smaller  covariance  at  shorter  time  periods  probably 
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comes  from  the  fact  that  stronger  winds  cause  more  mixing 
and  thus  decrease  AT  and  Ag. 

Results  from  NFS  cruise  data  show  the  SAM  (using 
average  wind  speeds,  temperatures,  relative  humidities  and 
pressures  for  the  entire  period  of  the  exercise)  gave  heat 
^®si-ilts  that  are  ~5%  (6%  LHF  and  3%  SHF)  larger  than 
when  using  the  sampling  method  (Table  5.8). 


Exercise 

Region 

LHF 

SHF 

THF 

CEAREX 

Greenland  Sea 

1.095 

1.000 

1.037 

SIZEX 

Greenland  Sea 

1.054 

1.071 

1.062 

NORCSEX  88 

Norwegian  Sea 

1.052 

1.077 

1.062 

NORCSEX  91 

Norwegian  Sea 

1.068 

0.997 

1.055 

CEAREX 

Barents  Sea 

1.130 

0.990 

1.086 

SIZEX 

Barents  Sea 

0.941 

1.005 

0.961 

Table  5.8.  Ratio  of  heat  flux  (W/m^)  calculated  using  time 
averages  of  meteorological  variables  versus  using  the 
sampling  method. 
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2 .  Stability  Dependence 

If  stability  is  neglected,  then  one  must  assume  neutral 
or  near  neutral  conditions  and  buoyant  convection  would  be 
negligible.  However,  neutral  stability  is  rare  and  unstable 
conditions  prevailed  during  the  NPS  exercises.  The  effect 
of  stability  is  applied  through  use  of  the  stability 
functions  and  'PJ  defined  in  Equations  4.11  and  4.12  for 
unstable  stratification  and  Equations  4.14  and  4.15  for 
stable  stratification.  For  unstable  stratifications 
(negative  z/L)  both  stability  functions  ('P)  are  complex 
functions  of  stability,  and  their  values  will  be  relatively 
small  and  negative  thus  yielding  larger  LHF  and  SHF  than 
under  neutral  conditions.  For  stable  stratification 
(positive  z/L)  ,  ^  is  a  linear  function  of  stability  ('P  =  - 
5z/L)  and  will  always  yield  smaller  LHF  and  SHF  than  neutral 
conditions . 

The  method  used  by  Bunker  accounts  for  stability.  HC 
however,  do  not  explicitly  include  stability  in  their 
calculations . 

3 .  Transfer  Exchange  Coefficients 

The  average  Cg  and  Cj,  for  the  NPS  Barents  Sea  exercises 
were  1.35x10-^  and  1.13x10-^  ,  respectively.  If  Bunker's 
nomogram  for  transfer  coefficients  were  used  instead,  that 
value  would  be  1.58xl0'\  This  17%  and  40%  increase  in  the 
exchange  coefficients  would  lead  to  a  corresponding  increase 
in  LHF  and  SHF.  Bunker’s  and  C^n  values  for  near-neutral 
conditions  were  strongly  dependent  upon  wind  speed  and  AT. 
This  dependence  on  wind  speed  is  however  unwarranted  as 
there  does  not  seem  to  be  conclusive  experimental  evidence 
for  variation  of  Cg  and  with  wind  speed  (Smith,  1988) , 

HC  heat  flux  errors  would  result  from  the  use  of 
constant  exchange  coefficients,  vice  coefficients  that  vary 
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with  wind  speed,  temperature  gradient,  humidity  gradient  and 
stability.  HC  fixed  and  at  1.5x10'^  for  all  regions. 

In  contrast,  the  coefficients  used  here  varied  by  region 
with  Cj,  having  a  range  of  1.0x10'^  to  1.25x10'^  with  an 
average  of  1.13x10“^  and  Cg  having  a  range  of  1.2x10"^  to 
1.5x10"'’  with  an  average  of  1.35x10"^.  If  calculating  fluxes 
for  our  exercises,  this  variation  would  cause  HC  to 
overestimate  SHF  nearly  45%  in  the  Barents  Sea  and  27%  in 
the  Norwegian  and  Greenland  Seas  and  to  overestimate  LHF  25% 
in  the  Barents  Sea  and  9%  in  the  Norwegian  and  Greenland 
Seas . 


4 .  Other  Factors 

As  discussed  in  Chapter  II,  the  VT  data  are  divided 
into  five  degree  latitudinal  sections  and  no  distinctions 
are  made  for  longitude.  Based  on  this  I  have  labeled  the  VT 
75N  -  SON  region  as  the  Greenland  Sea,  their  75N  -  70N  as 
the  Barents  Sea  and  their  70N  -  65N  as  the  Norwegian  Sea. 
These  distinctions  are  not  absolutely  accurate  in  a 
geographic  sense,  but  they  are  fitting  with  the 
climatological  stations  of  Gorshkov  and  Bunker  and  with  our 
ship  cruise  data. 

Haikkinen  and  Cavalieri  used  gridded  data  for  the 
calculation  of  heat  flux.  Thus,  comparison  of  locations 
should  be  easily  accomplished.  However,  their  data  was  for 
only  one  year  (1979),  which  will  likely  cause  errors  when 
comparing  with  different  years.  Other  errors  might  occur  by 
use  of  model  data  vice  observed.  If  the  model  was  biased 
toward  higher  winds  or  colder  air  temperatures  then  larger 
heat  fluxes  would  be  shown.  Additionally,  the  coarse  model 
grid  (250km)  does  not  have  the  resolution  to  accurately 
portray  the  large  horizontal  gradients  of  temperature  and 
wind  speed  that  occurs  at  the  ice  edge  (Personal 
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communication  Guest,  1995),  This  may  in  part  account  for 
their  large  flux  gradients  near  the  ice  edge. 

D.  POWER  SPECTRAL  DENSITY  OF  HEAT  FLUXES 

Power  spectral  densities  (psd)  of  heat  flux  for  all 
exercises  were  also  examined.  The  power  spectral  density 
function  estimates  the  signal  of  heat  flux  as  a  function  of 
frequency.  Windowing  was  performed  using  a  Hanning  window. 
The  sampling  frequency  that  varied  for  each  exercise  based 
on  the  record  length  was  used  for  scaling  plots.  One  half 
of  the  sampling  frequency  (1  /  600  sec  )  would  be  the 

highest  resolvable  frequency  (Nyquist  frequency) .  The 
lowest  resolvable  frequency  would  be  1  /  the  total  time  of 
each  exercise. 

All  the  power  spectral  density  plots  (Figures  5.13  - 
5.18)  exhibited  a  peak  in  power  between  1x10'^  -  5x10'^  Hz, 

corresponding  to  1-2  days.  This  illustrates  that  most  of 
the  power  is  at  the  synoptic  period.  The  typical  slope  of 
the  psd  function  at  the  higher  frequencies  (above  10'^  Hz) 
corresponding  to  periods  less  than  one  day  was  -1.5  in  log- 
log  space. 

E.  ACCURACY  OF  HEAT  FLUX  MEASUREMENTS 

The  largest  source  of  uncertainty  for  SHF  and  LHF 
measurements  is  sensor  accuracy  (Blanc  1987) .  Using  the 
manufacturers  published  sensor  errors  along  with  field 
experience  (Table  5.9)  an  error  estimate  (Tables  5.10  and 
5.11)  caused  by  sensor  inaccuracies,  ship  contamination  and 
^sdiation  contamination  for  our  exercises  was  determined  in 
a  root  mean  square  manner  (Equations  5.1  and  5.2). 
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Air  temperature  Error 


1.0  C 


Sea  temperature  Error 

1.0  C 

ASTD  Error 

1.4  C 

Relative  Humidity  Error 

5% 

Abs .  Humidity  Error  at  -10  C 

0.1  g/Kg 

Abs .  Htimidity  Error  at  0  C 

0.2  g/Kg 

Abs.  Humidity  Error  at  10  C 

0.4  g/Kg 

Wind  Speed  Error 

0.3  m/  s 

Table  5.9  Possible  measurement  errors  of  data  used  in  heat 
flux  calculations. 


(p  CpC’*  (  C/)  )  2+  (p  CpCj,  (AD  )  2  (5.1) 

=v'(pl.^Cg(Ag„^)  (JJ)  )2+  (pL^,C^(Ag)  (5.2) 
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5  m\s 

10  m\s 

15  m\s 

20  m\s 

25  m\s 

0  c 

10.3 

20.6 

30.8 

41.1 

51.4 

-5  C 

10.5 

20.7 

30.9 

41.2 

51.5 

-10  C 

11.4 

21.2 

31.2 

41.4 

51.7 

-15  C 

12.2 

21.6 

31.5 

41.9 

51.8 

-20  C 

13.5 

22.4 

32.0 

42.0 

52.2 

Table  5.10  Possible  error  in  sensible  heat  flux  (W/m^) 
measurement  as  a  function  of  wind  speed  and  ASTD. 


5  m\s 

10  m\s 

15  m\s 

2  0  m\s 

25  m\s 

0  g\Kg 

8.8 

17.6 

26.3 

35.1 

43.9 

-1  g\Kg 

8.9 

17.7 

26.3 

35.1 

43.9 

-2  g\Kg 

9.2 

17.8 

26.4 

35.2 

44.0 

-3  g\Kg 

9.6 

18.0 

26.6 

35.3 

44.1 

Table  5.11  Possible  error  in  latent  heat  flux  (W/m^) 
measurement  as  a  function  of  wind  speed  and  air-sea  humidity 
dif erence . 


The  random  SHF  and  LHF  errors  for  average  values  of 
wind  speed  (10  m/s),  ASTD  (-5  C)  and  air-sea  humidity 
difference  (-1.5  g/Kg)  from  the  NPS  data  set  would  be  21 
W/m^  and  18  W/m^  respectively.  Therefore  the  random  THF 
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error  would  be  about  27  W/m^,  when  added  in  a  root  mean 
square  fashion.  Random  errors  can  however  be  reduced  by 
increasing  the  number  of  averaging  operations  done.  This 
can  be  accomplished  without  changing  method  by  taking  more 
measurements.  Systematic  (non-random)  errors  are  also 
present  and  cannot  be  reduced  by  averaging.  The  worst  case 
non-random  error  for  the  NFS  exercises  would  be  39  W/m^. 
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Figure  5.2  Histogram  of  total  turbulent  heat  flux  (W/m*)  for 
SIZEX  in  the  Greenland  Sea,  January  1992. 
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Julian  date 


Figure  5.3  Time  series  of  sensible,  latent  and  total 
turbulent  heat  flux  (W/m‘)  for  CEAREX  in  the  Greenland  Sea 
March  1989. 


#  of  occurrences 
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Figure  5.4  Histogram  total  turbulent  heat  flux 
CEAREX  in  the  Greenland  Sea,  March  1989. 


(W/m") 


for 


600 
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julian  date 


Figure  5.5  Time  series  of  sensible,  latent  and  total 
turbulent  heat  flux  (W/m^)  for  NORCSEX88  in  the  Norwegian 
Sea,  March  1988. 


#  of  occurrences 


121 


Figure  5.6  Histogram  of  total  turbulent  heat  flux  (W/m^)  for 
NORCSEX  in  the  Norwegian  Sea,  March  1988. 
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Figure  5.8  Histogram  of  total  turbulent  heat  flux  (W/m‘)  for 
NORCSEX  in  the  Norv/egian  Sea,  November  1991. 
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Figure  5.9  Ti; 
turbulent  hea 
F  ebruarv-Mar c 


Figure  5.10  Histogram  of  total  turbulent  heat  flux  (W/m") 
for  CEAREX  in  the  Barents  Sea,  March  1989. 
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Figure  5.11  Time  Series  of  sensible,  latent  ant  total 

turbulent  neat  flux  (W/m‘)  for  SIZEX  in  the  Barents  Sea 
March  1992.  dj-cnui,  oea, 


Figure  5.12  Histogram  of  total  turbulent  heat  flux  (W/m") 
for  SIZEX  in  the  Barents  Sea,  March  1992. 
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frequency 


Figure  5.16  Power  spectral  density  of  sensible 
line),  latent  (dashed  line)  and  total  turbulent 
heat  flux  (W/m"')  for  N0RCSEX91  in  the  Norwegian 
November  1991. 


(dotted 
(solid  1 
Sea, 


frequency 


Figure  5.17  Power  spectral  density  of  sensible  (dotted 
line),  latent  Jdashed  line)  and  total  turbulent  (solid  line) 
heat  flux  (W/iri-)  for  CEAREX  in  the  Barents  Sea,  Februarv- 
March  1989. 


VI .  CONCLUSIONS 


This  thesis  examined  the  exchange  of  heat  between  the 
ocean  and  atmosphete  in  the  Greenland,  Norwegian  and  Barents 
Seas .  The  bulk  method  was  used  to  calculate  the  turbulent 
heat  fluxes  in  the  winter  based  on  bulk  data  gathered  aboard 
ships.  The  results  showed  the  largest  positive  fluxes 
occurring  in  the  Greenland  Sea,  followed  by  the  Norwegian 
Sea  and  the  smallest  fluxes  in  the  Barents  Sea.  The 
Greenland  Sea  is  also  a  location  of  strong  convective 
overturning  and  deep  water  formation  in  the  Arctic. 

The  values  of  the  turbulent  heat  flux  were  generally 
found  to  be  consistent  with  previous  climatological  results. 
However,  HC  flux  estimates  in  proximity  to  the  ice  edge  were 
much  larger  than  those  observed  during  the  NFS  cruises . 
Specifically,  they  estimated  mean  turbulent  heat  fluxes  well 
in  excess  of  400  W/m^  and  even  as  high  as  600  to  700  W/m^  in 
portions  of  the  Greenland  and  Barents  Seas .  Based  on  the 
NFS  data  and  in  comparison  to  the  climatological  works  in 
the  regions,  the  estimates  of  HC  are  extremely  high.  In 
fact  their  mean  values  for  some  regions  of  the  Greenland  and 
Barents  Seas  are  nearly  equivalent  to  maximum  fluxes  found 
in  this  study.  Differences  in  their  heat  flux  values  and 
those  in  this  thesis  may  be  due  in  part  to  larger  winds  in 
the  HC  data.  Also,  HC  used  moisture  and  heat  exchange 
coefficients  about  11%  and  33%  larger  than  our  measured 
f i ^i ®^^t s ,  explaining  some  of  the  differences. 

A  large  uncertainty  in  the  flux  measurements  of  Bunker 
and  HC  resulted  from  the  use  of  questionable  heat  and 
moisture  exchange  coefficients.  Bunker  used  a  stability  and 
wind  dependent  exchange  coefficients  that  appeared  to  give 
values  about  40%  larger  than  for  the  method  of  Smith  (1988). 
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HC  used  constant  exchange  coefficients  that  would  be  11-33% 
larger  than  those  of  Smith  (1988)  for  the  NPS  data  sets  used 
in  this  thesis. 

Finally,  the  averaged  values  of  the  bulk  meteorological 
parameters  can  be  used  in  the  bulk  formula  to  obtain 
averaged  sensible  and  latent  heat  fluxes  within  5%  of  the 
actual  values.  These  results  validate  the  classical  monthly 
averaged  method  (SAM)  employed  by  VT,  Gorshkov  and  others  in 
the  Arctic  Seas.  However,  other  errors  in  their  method  may 
cause  inaccurate  results.  A  principal  advantage  of  the 
classical  method  is  speed  of  processing.  After  the  raw  data 
has  been  quality  controlled,  an  average  heat  flux  can  be 
computed  in  seconds  via  the  classical  method  compared  with 
one  hour  or  more  with  the  sampling  method. 
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APPENDIX.  [MATLAB  PROGRAM  FOR  COMPUTING  HEAT  FLUX] 


matlab  4.0  program  run.m 
updated  Feb  1995 

based  on  smith  (jgr,v  93, nO  cl2,  15467-15472,  1988] 
calculate  the  following  fluxes 


shf 
Ihf 
hf 

based  on 

rho  atm  denstiy  based  on  sfc  temp 

ustar  see  below 

tstar  see  below 

l_v  specific  heat  set  constant  at  1004  j  kg-1  k-1 

c_p  latent  heat  of  vaporiz  const  at  2.50e+6  j  kg-1 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%^^^^^^^^^^^^^^^ 
zu=15.;  % 

zt=15.;  % 

load  norx88; 
myf ile=norx88 ; 

[n  m] =size (myf ile) ; 
for  j=l:n; 

min(j)  =  (j-1)  *  10; 
date ( j ) =myfile { j , 1)  ; 
time  ( j  )  =iTyfile  ( j  ,  2 )  ; 
utrue ( j ) =myfile ( j ,7); 
tair(j)=iryfile(j,3)  ; 
tsfc ( j ) =myfile ( j , 4)  ; 
rh ( j ) =myf ile ( j , 5) ; 
p ( j ) =myfile ( j , 6) ; 
wdt ( j ) =myf ile ( j , 8) ; 
lat ( j ) =myfile ( j , 9) ; 

Ion ( j ) =myfile ( j  ,  10)  ; 

jd(j)  =  67  +  (1200  /  1440)  +  (  min(j)  /  1440); 

wind ( j ) =utrue ( j ) ; 

astd(j ) =tair ( j ) -tsfc ( j ) ; 

press { j ) =p ( j ) ; 

relhum ( j ) =rh { j ) ; 

sst ( j ) =tsfc (j); 

temp ( j ) =tair ( j ) ; 

utrue=utrue ( j ) ; 

tair=tair ( j )  ; 

tsf c=tsfc ( j ) ; 

rh=rh( j )  ; 

P=P  ( j  )  ; 

i  calculate  bulk  ustar, tstar , gstar 


sensible  heat  flux 
latent  heat  flux 
(total)  heat  flux 


/  1440 


min(j)  /  1440; 


input 


utrue 

tair 

rh 


wind  speed  (m/s)  at  zu 
air  temperature  (centigrade) 
relative  humidity  (%)  at  zt 
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tsfc  surface  temperature 

zu  measurement  level  of 

zt  measurement  level  of 

p  atmospheric  pressure 


centigrade) 
utrue  (centigrade] 
tair  and  rh  (m) 
at  zt  (mb) 


output 


ustar  friction  velocity  (m/s) 
tstar  scaling  temperature  (k  or  c) 
gstar  scaling  specific  humidity  (g/kg)  not  (g/g) ! 

1  monin-obukov  length  scale  (m) 

calculates  tstar  based  on  a  constant  chn,cen 
zo  base  on  charnocks  relation 

this  follows  smith (1988)  as  close  as  psby  except  different 
virtual  temperature  calculation  is  used  (smith  was  wrong) 
the  surface  is  assumed  to  be  saturated, 
set  constants 

zl0=10;  %!  reference  height  (not  measuremement  height) 

next  two  are  based  on  value  at  zlO 


ctn=l . Oe-3 ; 
cen=l .2e-3 ; 
k= .  4  ; 
g=9.8; 
a=0.011; 
b=0.11; 
nu=l . 4e-5 ; 
gamma=0 . 0096 , 


heat  and  buoyancy  flux  xfer  param  at  zlO 
1  humidity  flux  transfer  parameter 
I  von  karmen's  constant 
.  gravity 

charnocks ' s  constant  from  smith (1980) 
smooth  flow  constant  from  businger ( 1973 ) 
dynamic  viscosity  of  air 
adiabatic  lapse  rate 


^%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%% 
calculate  parameters  not  requiring  iteration 
calculate  saturation  mixing  ratio,  qsat  (g/kg) , 
based  on  temperature  (c)  and  press  (mb)  from  stull  (1988) 
it  is  only  valid  near  the  surf  ace.  use  loew's  polynomial  for 
more  accurate  or  upper  level  estimates  of  qsat. 
psfc=p+0.116  *  zt;  %  sfc  press  (mb)  based  on  standard  atms 
es=6 . 1078*exp (17 . 2694  *tsfc/  (tsfc+237 .3 ) ) ; %  vap  p (mb)  tsfc 
qsat=622.0  *  es/  (psfc-es);  %  (g/kg) 

esa=6.1078*exp(17.2694  *tair/  ( tair+237 . 3 ) ) ; %vap  p (mb)  tsfc 
qair=622.0  *  esa/  (p-esa) ;  %  (g/kg) 

q=qair  *  rh/100;  %  mixing  ratio  at  zt  (g/kg) 
qsfc=qsat;  %!  mixing  ratio  at  surface  assumed  staturated 
theta=tair+273 .16+gamma  *  zt;  %!  potent  temp  (k)  at  zt 
the  following  is  based  on  stull (1988) ;smith(1988)  is  wrong 


thetav=theta* (1 . 0+0 . 61e-3*q)  ;  %!  virt  potent 
tvsfc= (tsfc+273 . 16) * (1 . 0+0 . 61e-3*qsfc) ;%  vi] 
)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


:!  virt  potent  temp  (k)  at  zt 
-3*qsfc);%  virt  pot  temp  sfc 


initialization  (these  values  will  change) 
count =0 ; 
l=1.0el0; 
zo=l . Oe-4 ; 
ustar=utrue  *  .036; 

find  ustar,  tstar  loop;  exit  after  convergence  or  40  cycles 
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%  count  is  var  for  the  iterations 
if  utrue  >  0.2; 

%  calculate  an  init  ustar 
ustarl=ustar ; 
zul  =  zu  /  1; 
ztl  =  zt  /  1; 
ratio  =  zu  /  zo; 

calculate  the  integral  diabatic  heating  term 
for  momentum,  we  have  used  values  from  dyer (1974) 
unstable 
if  zu/1  <  0 . ; 

x=(l.-16.*(zu/l)  )''.25; 

!im=2  .  *logm(  (1  .+x)  /2  . )  +logm(  (l.+x*x)  /2 )  -2*atan(x)  +3.14159/2  ; 
stable 
else 

psim  =  -5  *  (zu/  1) ; 
end 

for  temperature,  based  on  large  and  pond  (jpo  v.  12,1982) 
unstable 

if  zu/1  <  0.; 

x=(l.-16.*(zu/l))'^.25; 
psit=2 . *logm( (l.+x*x)/2) ; 

%  stable 
else 

psit  =  -5  *  (zu/  1  ) ; 
end 

if  ratio  <5.0  ;  %  prevent  underflow 
count =40 ; 
end 

zot=zl0*exp (-k*k  /  (ctn  *  logm  (zu/zo)));  % 

tstarv=  (thetav-tvsfc)  *  k  /  (logm(zt  /  zot)-psit); 
ustar=utrue*  k  /  (logm  (zu  /  zo)-psim); 
zc=a  *  ustar*  ustar/g;  %  charnocks  formula 

zs  -b*nu  /  ustar;  %  smooth  formula 

zo=zc+zs ; 

if  abs(tstarv)  <  l.Oe-20; 

l=1.0el0; 

else 

l=ustar*  ustar*  tvsfc/  (g*k  *  tstarv) ; 
end 

%  do  iterations  to  find  ustar  convergence 
if  abs ( (ustar-ustarl) /ustar)  >  .0005; 
for  i=l:40; 
ustarl=ustar ; 
zul  =  zu  /  1; 
ztl  =  zt  /  1; 
ratio  =  zu  /  zo; 

■6  this  function  calculates  the  integral  diabatic  heating  term 
■6  for  momentum,  we  have  used  values  from  dyer  (1974) 
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%  unstable 

if  zu/1  <  0 .  ; 

x=(l.-16.* (zu/1) ) ^.25; 

psiin=2  .  *logm(  (l.+x)  /2  .  )+logm(  (l.+x*x)  /2 )  -2*atan (x) +3 . 14159/2 ; 
%  stable 

else 

psim  =  -5  *  (zu/  1) ; 
end 

for  temperature,  based  on  large  and  pond  (jpo  v.  12,1982) 
unstable 

if  zu/1  <  0 . ; 

x=(l.-16.*(zu/l) )".25; 
psit=2 . *logm( ( 1 . +x*x) /2 ) ; 

%  stable 

else 

psit  =  -5  *  (zu/  1  ) ; 
end 

if  ratio  <5.0  ;  %  prevent  underflow 
count =40 ; 

end 

zot=zl0*exp {-k*k  /  (ctn  *  logm  (zu/zo)));  % 

tstarv=  (thetav-tvsfc)  *  k  /  (logm(zt  /  zot)-psit); 
ustar=utrue*  k  /  (logm  (zu  /  zo)-psim); 
zc=a  *  ustar*  ustar/g;  %  charnocks  formula 

zs  =b*nu  /  ustar;  %  smooth  formula 

zo=zc+zs ; 

if  abs(tstarv)  <  l.Oe-20; 

l=1.0el0; 

else 

l=ustar*  ustar*  tvsfc/  (g*k  *  tstarv) ; 
end 

end  %  is  for  i  at  40 
end  %  for  abs  (  )  > 

end  %  when  u< . 2 

o  post  loop  calculations 

if  count  <  40  &  abs ( (ustar-ustarl) /ustar)  <  .0005; 
tstar= ( tair-tsfc)  *  k  /  (logm(zt/  zot)-psit); 
gstar= (q-qsfc)  *k  /  (logm(zt/  zot)-psit); 
else  %  too  stable  loop  does  not  converge 
ustar=0 . ; 
tstar=0 . ; 
qstar=0 . ; 

1  =  0.  ; 
end 

rho=psfc/ (2 . 87*tvsfc)  ;  %  calc  var  density  with  ideal  gas  law 
shf ( j ) =-rho*1004  *  ustar*  tstar; 

Ihf ( j ) =-rho*2 .5e6  *  ustar*  qstar/1000 ; 1000 
hf ( j ) =shf ( j ) +lhf ( j ) ; 
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cd  ( j  )  =  (uustar  ( j  )  /  utrue)  2; 
cdn{j)=  (1  /  sqrt  (cd(  j  )  ) +psiin/ .4)  (-2); 

fprintf ( 'working  on  observation  number  %g\n',j) 
end  %%  stop  counting  j  here 
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